Das Ondas Contínuas para Matrizes Discretas
A vibração de uma corda ou membrana é governada pela equação de onda:
$$\rho(x) \frac{\partial^2 v}{\partial t^2}(x, t) = \frac{\partial}{\partial x} \left[ p(x) \frac{\partial v}{\partial x}(x, t) \right]$$
Para encontrar a solução $v(x, t) = \sum_{k=0}^{\infty} c_k u_k(x) \cos \sqrt{\lambda_k}(t - t_0)$, devemos resolver o sistema de Sturm-Liouville:
$$\frac{d}{dx} \left[ p(x) \frac{du_k}{dx}(x) \right] + \lambda_k \rho(x) u_k(x) = 0$$
A discretização do operador leva a equações matriciais como $Aw = -0.04 \frac{\rho}{p} \lambda w$. Para uma matriz tridiagonal $4 \times 4$, $p(\lambda)$ é gerenciável. No entanto, à medida que a malha se refina ($n$ aumenta), enfrentamos dois obstáculos:
- Limite de Abel-Ruffini: Não existe solução algébrica para as raízes de polinômios onde $n \ge 5$.
- Sensibilidade ao Arredondamento: Em sistemas de alta dimensão, uma mudança na casa decimal $10^{-10}$ de um elemento pode deslocar os autovalores por ordens de grandeza (fenômeno do polinômio de Wilkinson).
Necessidade Numérica e Bibliotecas Profissionais
Bibliotecas numéricas profissionais (IMSL, NAG) evitam polinômios característicos brutos. Em vez disso, utilizam rotinas iterativas para aproximação:
- Biblioteca IMSL: Utiliza mínimos quadrados lineares, splines cúbicas e transformadas rápidas de Fourier.
- Biblioteca NAG: Emprega aproximação polinomial por mínimos quadrados e ajustes no sentido de $l_1/l_{\infty}$.
Ao aproximar autovalores para o sistema $\lambda_i = 1 + 4\alpha\left(\sin \frac{\pi i}{2m}\right)^2$, dependemos de mínimos quadrados discretos e descoberta iterativa em vez de busca de raízes.